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Abstract. - We theoretically study divergent fluctuations of dynamical events at non-ergodic 
transitions. We first focus on the finding that a non-ergodic transition can be described as a saddle 
connection bifurcation of an order parameter for a time correlation function. Then, following the 
basic idea of Ginzburg-Landau theory for critical phenomena, we construct a phenomenological 
framework with which we can determine the critical statistical properties at saddle connection 
bifurcation points. Employing this framework, we analyze a model by considering the fluctuations 
of an instanton in space-time configurations of the order parameter. We then obtain the expo- 
nents characterizing the divergences of the length scale, the time scale and the amplitude of the 
fiuctuations of the order parameter at the saddle connection bifurcation. The results are to be 
compared with those of previous studies of the four-point dynamic susceptibility at non-ergodic 
transitions in glassy systems. 
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Slow relaxation behavior is observed generally in glassy 
Systems, such as super-cooled liquids, dense colloidal sus- 
pensions, foams, emulsions and granular materials. In 
particular, it has been found that for some systems, the 
time correlation function exhibits a plateau over a cer- 
tain time interval. Although the question of whether or 
not the plateau becomes infinitely long at some transition 
temperature has not been definitely answered, there is ex- 
perimental evidence that the plateau does become very 
long as a particular temperature is approached. This dis- 
tinctive phenomenon is called a non-ergodic transition [1] . 
Recently, in addition to studies of the singular behavior 
of glassy systems, there have been investigations of the 
growth of the characteristic length scale and the enhance- 
ipent of the four-point dynamic susceptibility with the goal 
of elucidating the nature of the cooperative fluctuations of 
dynamical events responsible for the slow relaxation [2] . 

As one means for quantifying dynamical events in glassy 
systems, the following fluctuating quantity has been con- 
sidered [3,4]: 

qa{x,t) = 5p{x,t) d^ywa{x~y)Sp{y,0), (1) 



where Sp represents the deviation of the local density from 
its ensemble average and Wair) = cxp(— |rp/(2a^)) is an 
overlap function with size a. If qa{x,t) differs signifi- 
cantly from qa{x, 0), it can be concluded that a dynamical 
event leading to the de-correlation of the local density oc- 
curs during the time interval [0,t]. The ensemble average 
{qa{x, t)) is the time-correlation function of the local den- 
sity field. For the system with the spatially translational 



symmetry, it does not depend on x. Using the quantity qa, 
fluctuations of dynamical events are naturally quantified 
by the function C4(r,i) = {qa{r,t)qa{0,t)) - {qa{0,t)f. 
Then, it is natural to define the amplitude of fluctuations 
as Xiit) = j d'^rC4{r,t). It has been conjectured that 
X4(i) has a peak at some time t — t^, ^ \T — T^^'^ , that 
Xiit*) diverges as \T — Td""*", and that the correlation 
length r* of C4(r,t,) diverges as |T — Tc\~'^ , where Tc 
is the temperature at which a non-ergodic transition oc- 
curs [4]. 

There are several theoretical approaches to this phe- 
nomenon. A reliable one might be mode coupling theory, 
because it has been recognized that the theory can pro- 
vide accurate calculation for some experimental results of 
glassy systems. Recently, within a framework of mode 
coupling theory, the exponents of the divergences are cal- 
culated [5,6]. In contrast, as an alternative approach, one 
seeks for another framework that is independent of spe- 
cific properties of glassy systems. For example, a space- 
time thermodynamic formalism has been proposed with 
this motivation [7]. In this Letter, we present a new the- 
ory in the latter approach. The framework we propose is 
analogous to the Ginzburg-Landau theory of critical phe- 
nomena [8]. Within this framework, we calculate the ex- 
ponents of the divergences by applying a singular pertur- 
bation method to a path integral expression. 

Framework:. — Let us begin by recalling that the ef- 
fective Hamiltonian in Ginzburg-Landau theory describes 
a pitchfork bifurcation in the simplest form; that is, the 
Ginzburg-Landau potential possesses a double well form 
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below a critical temperature. This simplicity is responsible 
for the universality of the results of the Ginzburg-Landau 
theory. With this property of Ginzburg-Landau theory 
in mind, the first problem we consider here is to identify 
the type of the bifurcation occurring at non-ergodic tran- 
sitions. In fact, this problem is solved in Ref. [9]: The 
behavior of the time correlation function is determined by 
an equation, dfcp = fe{4>, dt(j>) that exhibits a saddle con- 
nection bifurcation [10]. Here, the order parameter </> is 
defined as follows. First, the equation for the time corre- 
lation function C{r,t){= {Sp{r,t)5p{0,0))) is derived as 
dfC = !F{C) from a many-body Langevin equation, un- 
der the assumption that third- and fourth-order cumu- 
lants of density fluctuations are ignored in d^C. We then 
find that the linear operator 8T/8C around C — Q has 
a zero eigenvalue at some temperature T = Tq with de- 
creasing T. Using the zero eigenvector $00 (^) that satisfies 
8J- /5C ■ $00 = at r = To, wo define the order parameter 
as the amplitude of $oo('") for C{r,t). See Ref. [9] 
for the detail. 

Here, wc! note that the results presented in this work 
are independent of the form of /, as long as it exhibits 
a saddle connection bifurcation at e = with decreasing 
e, as wc find below. For this reason, instead of employ- 
ing the complicated form of fe{4>,dt(j)) derived for a spe- 
cific model of many-body Langevin systems in Ref. [9], 
we assume the simple form f^((j),dt4>) = —d,f,Uf_{<j)) with 
U^{(l)) = -^2 _ _^ ^ .^jjgj.g g > [11]. The exis- 
tence of the saddle connection bifurcation in this model 
is easily seen from the energy conservation condition 
{d(j)/dt)'^ /2 + U^{(j)) = for trajectories of (j){t) that satisfy 
a boundary condition (j){oo) = 0. This equation exhibits 
the following behavior. When e is small, (p first approaches 
0=1 from 0(0) = 00 > 1) and then after a delay, it be- 
gins to approach the origin, = 0. As e — > 0, the time 
interval during which remains near 1 diverges. This type 
of bifurcation is termed a saddle connection bifurcation, 
because the two fixed points 0=1 and = 0, which form 
saddles in the two-dimensional phase space (0, 9i0), are 
connected by a trajectory for the case e = 0. There is a 
clear qualitative correspondence between 0(t) [Q <t < 00) 
and a time correlation function exhibiting a non-ergodic 
transition. 

Now, we employ the basic idea of Ginzburg-Landau the- 
ory. First, in order to describe large-scale spatial vari- 
ation of the order parameter 0(t), we replace the order 
parameter 0(t) with an order parameter field (j>{x,t) [8]. 
When fluctuations of are ignored, as the simplest de- 
scription of spatially large-scale behavior of the order pa- 
rameter field, we assume that 4'{x,t) obeys the evolution 
equation dfcf) = fc{4>) + Dl^^cj) under the boundary con- 
dition 0(a;, 00) = 0, where DA'^cp represents the diffusive 
coupling associated with 0^0, and we write f^{(j),dt4') as 
/e(0), for simpHcity. Here, = (V • V)^. Note also 
that periodic boundary conditions are first assumed for 
the system with a size L, and then the limit L — > 00 is 
taken. Next, we take account of the fluctuation effects of 



(j){x, t) by considering space-time configurations of 0(a;, t) 
with (x e and < t < 00), which we denote by [0]. 
Then, let f ([0]) be the probability measure for [0]. Until 
now, we have not derived the accurate form of ^([0]) on 
the basis of microscopic models. Nevertheless, relying on 
a phenomenological argument, we can reasonably assume 
a form of ^'([0]) in the following manner. 

First, we divide time t as tk = kAt, with k = 
1,2- ••. Denoting a collection of space configurations of 
{(j){x,t),dt(j>{x,t)) at an instant time t as u(f), we assume 
that there is a time interval At such that Tm <C At <c tm, 
where Tm is the longest time scale of fast variables elimi- 
nated in deriving the order parameter equation from mi- 
croscopic models, and tm is the characteristic time scale 
of the time evolution of u. We then consider the transition 
probability T{u u') during a time interval At. Now, as 
in the case of the Onsager-Machlup theory for stochastic 
processes describing the time evolution of thermodynamic 
variables [12], the assumption mentioned above suggests 
the following two properties: (i) T{u u') can define a 
Markov process and (ii) T{u — » u') possesses the large 
deviation property with respect to time. From the first 
property (i), we can express P{[(j>]) as 



(2) 



P([0]) ~ Po{uo)T{uo Ui)T{ui -> U2) ■ 



in the discretized description, where Uk = u{kAt) and 
Pq{uq) is the stationary distribution function of u at an 
instant time. Further, from the second property (ii), T can 
be expressed as T{u u') ~ exp(— At/(M'|M)) with a large 
deviation functional I{u'\u). Here, I(u'\u) — should 
provide the discretized form of the deterministic equation 
9^0 = /e(0) + DA^cf). Thus, approximating I{u'\u) with 
the second order polynomial of u' , we assume 



X 



[(9*0)' - 9*0 - Ai(/,(0) + i)A20)]V(A^)^ (3) 

where B is a constant. From (2) and (3), taking the limit 
At/TM — > 0, we can express P([0]) as 



n[<^]) = |exp 



with 



ll^^dtj d'x{d^<p - fM - DA^f 



(4) 



(5) 



under the boundary condition 0(a;,oo) = 0. Here, in- 
stead of considering Po(mo)! we fix the initial condition 
0(a;,O) as 0o that is determined from the equal-time cor- 
relation function. Note that the so-called Jacobian term 
is not written in (5), because it is independent of [0] in 
the present problem [13]. 

Although the model (4) with (5) is not justified by the 
analysis of microscopic models, we expect that this model 



T^.0 



Divergent fluctuations at non-ergodic transitions 



can capture essential aspects of fluctuations near the sad- 
dle connection bifurcation. With this expectation, we cal- 
culate the quantity {^{x, t)) and 

t) = {<t>{x, t)(P{0, t)) - {<l>ix, t)) {<l>{0, t)) (6) 

employing (4), and then from this we compute X<p{t) = 
J d^xCtp{x,t). Because {(f){x,t)) determines the behavior 
of {qa{x,t)), we believe that (j>{x,t) characterizes dynam- 
ical events responsible for slow relaxation. We thus con- 
jecture; that X'pif) has a peak at some t = and that 
the system displays divergences of the forms t* ~ e~^, 
X<p{t*) — e"''') and r* ~ e~'^, where r* is the correlation 
length defined from Cci){x,tt). In the following, we demon- 
strate that indeed such divergences do occur and that we 
have the values ( = 2, j — 3/4 and v = 1/4. 

Analysis:. — The probability distribution (4) can be 
realized as the stationary distribution of the four dimen- 
sional field <j){x, t) that obeys a fictitious stochastic pro- 
cess. Introducing a fictitious time s, we can construct such 
a stochastic model as 



ds(j){x,t;s) = 



5F 



6(j){x, t; s) 



-^{x,t;s), 



(7) 



where 



{^{x,t;s)^{x' ,t'; s')) = 2BS^{x-x')S{t-t')S{s- s'). (8) 

It is easily checked that the s-stationary distribution func- 
tion of <j>{x,t) for this model is exactly expressed by (4). 
Then, using the operator Lg{(j)) = df — d^fe{(j)), we write 
(7) explicitly as 

ds<l> = -(L, - DA^)(dfct> - fM) - DA^<j>)+^. (9) 

Now, let (/)ci(^;e) be the solution of d^cj) = ft{4') with 
the boundary conditions (?!)ci(0; e) = 0o and (?ici(oo; e) = 0. 
Then, if i? = 0, all solutions of (9) approach this solution 
in the limit s oo. Furthermore, let and 0b (i) be 

special solutions of d^cf) = fo{<P) satisfying the conditions 
0,(— oo) = 1, (/)*(oo) = 0, (f)B{0) = 4'a, and cf)B{oo) = 1. 
Such solutions exist generally at saddle connection bifur- 
cation points. In the model under investigation, we have 
(/>*(t) = (1 — tanh(f/\/2))/2, which represents a kink solu- 
tion in the t direction, a kind of instanton [14] . Note also 
that (pnit) — Se~^*^ + 1 for Xt^ 1 with the constants A and 
S determined by the equation = /o(<^). Thus, for small 
e, the solution 4>ci{t', e) is close to (j)^{t — 9o) + (j)B{t) — l with 
a constant ^o- Based on these preliminary considerations, 
we find that, when e and B are small, it is convenient to 
express solutions of (9) as 

<p{x,t;s)=^(j)^{t-0{x,s)) + (l)B{t)-l + 'p{x,t;s), (10) 

where 0{x, s) represents the kink position, which depends 
on {x,s), and is assumed to be small when 6X^1 
and e <C 1. We then find that perturbations of the instan- 
ton decay exponentially, except for the Goldstone mode. 



corresponding to the translation of the instanton in the t 
direction [15]. It is known that in such a situation a sin- 
gular perturbation method is useful to derive an evolution 
equation for the kink position 9{x, s) under the infiuence 
of perturbations taking the forms of finite epsilon effects, 
(j)B{t) — 1 (see (10)), and noise [16]. 

This perturbative calculation is carried out by substi- 
tuting (10) into (9) and extracting the important terms 
under the conditions e 1, B <^ 1, 9X ^ 1 and that 
spatial variation of 6 is small. In order to make these as- 
sumptions explicit, we replace (pB — 1, e, V, and ^ with 
^{(jjB — 1), /^e, n^^^V and /x^^ and consider the expansions 

dgO = fifli + iJ^Q2 H and ip = iiipi + Ii^(p2 H , where 

/z is a small parameter, and fli and yi, are functions of 
the field 9{x, s). These replacements are determined in an 
essentially unique manner by the requirement that a sys- 
tematic perturbative expansion can be carried out. Then, 
collecting the terms proportional to /i, we obtain 

r2i5f(/), - Ll^i<pB - 1) = Lq4Lo*Vi - e/p(0*) - DA'^cj)^], 

(11) 

where Lq* — Lq{(I)^) and /p(0) = df^{(l))/de\^=Q. Because 
Lodt4>* = 0, the solvability condition for (11) leads to 
fii = [17]. Employing this condition, we obtain ipi. In 
a similar manner, the solvability condition for the equa- 
tion consisting of terms proportional to /x^ yields fl2, from 
which we derive 

BO 

ci — = ecae"^^ - e^cg - eCiDA'^e - ciD^A'^0 + S, (12) 

OS 

where wc have ignored terms non-linear in the gradient of 
0, such as (A6')^ and e~^^A9, and S satisfies 

{E{x, s)E{x', s')) = 2ciB5^{x - x')5{s - s'). (13) 

Note that the coefficients Cj are positive constants inde- 
pendent of e and can be calculated from (9). The first 
term on the right-hand side of (12) represents the interac- 
tion with the exponential tail of the solution 0b (i)- The 
second term in (12) determines the velocity of the steady 
propagation of the kink in a system without a boundary 
in the case e 0. 

Now, we investigate the s-stationary statistical proper- 
ties of 9 determined by (12) with (13). Wc first consider 
the case £) = 0, in which the field 9 behaves independently 
at each x. Thus, according to the central limit theorem, 
the spatial average of 9 over a domain C obeys a 
Gaussian distribution with average 0* and variance xe/l^l 
if the volume of the domain |0| is sufficiently large. In this 
case, the leading terms of f?* and xe in the limit e ^ can 
be calculated in an elementary manner and we obtain 



Xe 



B 

cse^ 
B^ 



(14) 
(15) 
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Next, we study the case in which £) is a small positive 
number. Specifically, defining 



ei 



d''xe"'''{0{x,s)-e^), 



we consider the s-stationary distribution function 

1 vail}) 



(16) 



(17) 



Although it is difficult to derive the exact form of y, we 
conjecture that 0(fe) with small |fc| obeys a Gaussian dis- 
tribution. Then, the fluctuation intensity of B(h = 0) 
should be equal to xo- Also, considering (12), we find 
that the length scale of d is simply determined by the 
balance among the second, third and fourth terms in the 
right-hand side when we focus on the regime :» 1. This 



leads to the scaling relation |fe| 



On the basis of these 



two considerations, we assume the simple form 



(27r)3 



2X0 



\0i{kr + ^e'k'\e,{k)\ 



(18) 

where D4 is a constant. Strictly speaking, (18) is valid 
only for small |fc|, but wc believe that it provides a good 
description of the large scale behavior in which we are 
interested. We therefore use this expression in order to 
evaluate {(p{x,t)) and {(p{x,t)^{x' ,t')) through the rela- 
tion (10). 

We calculate these quantities by first expressing 0*(t) 

as 

'^^ ■ (19) 



27r' 



Then, using {(p{x,t)) ~ {(f)^{t— 6{x))), we write 

{<l>{x, t)) J ^^4^)e«(*-^*) ^e— ^^(-)) . (20) 



Prom the Gaussian nature of the distribution function (17) 
with (18), we obtain 



-izOi {x) 



(21) 



with T,, 



cse 



-13/8 



, where C5 is independent of e. Thus, 



{4>{x, t)) can be expressed as a scaling form g{{t — 9*)/T^). 
Also, it is seen that (0(a;, t)) possesses a kink in the t di- 
rection whose average position and width are ^* ~ 0(e~^) 
and Tw — 0(e~^^/^), respectively. 

Next, we calculate {(j){x,t)4>{x' ,t)) by approximating it 
with {(j)^{t-e{x))(j)*(t-e{x'))). Then, using (19), we 
write 



{<j>{x,t)<i){x' ,t)) 



^exp 



i j d'^x"J\x")e{x") 



,(22) 



with J\x") = -zS^{x" -x)- z'6^{x" - x'). Performing 
the Gaussian integration of 9, we obtain 



^iz(t-0,)+iz'{t-9,) /g-iz'ei(x') 



exp 



-zz'TlG{ce\x-xy/^) 



(23) 



where cg is independent of e and G{r) = e^'' sin(r)/r. 
Specifically, in the long distance regime, for which we have 
|G'(c6|a; — a;'|e^/'^)| <^ 1, we can derive the approximate 
expression 

C4x,t) ~ [dt {<PiO,t))fGice\x\e'/^). (24) 

From this expression (23), we immediately find that the 
correlation length is proportional to e~^/^. Thus, we ob- 
tain ly = 1/4. Next, noting that the singular part of Xct>{t) 
originates from the long-distance part of C^{x, t), we can 
use the asymptotic form (24) in order to evaluate the ex- 
ponents C and 7. First, it is found that X(t>{i) has a peak 
at t = and thus, from (14), we obtain C, = 2. Fur- 
thermore, from the relation T^dtg{{t — 9*)/t^) ~ 0(e°), 
we obtain x^l^*) — e^'^^^, and hence 7 = 3/4. Finally, 
if the asymptotic formula (24) is valid on the correlation 
length scale, it suggests that the critical exponent 77 in the 
standard notation of critical phenomena is equal to —1. 

Conclusion:. — We have presented a phenomenolog- 
ical framework for determining the statistical properties 
of dynamical events at non-crgodic transitions. Within 
this framework, wc have obtained the results Q = 2, 
7 = 3/4 and u = 1/4, which are to be compared with 
the values obtained in numerical experiments (7/C — 0.4 
and i^/C — 0.22) [18] and the values calculated using the 
mode coupling theory (C ~ 2.3 for a Lennard-Jones fluid, 
7 = 1 [5,6], and u = 1/2 [5] and u = 1/4 [6].) [19]. It 
might be rather surprising that our theoretical values are 
close to the results of the mode coupling theory despite 
the great difference between the two frameworks. How- 
ever, because both theories are still at the mean field level 
in the sense that critical fiuctuations are not treated pre- 
cisely, it is not time to compare the results in detail. Sim- 
ilarly, with regarding the discrepancy with the results by 
numerical simulations, we wish to be cautious of a hasty 
judgment. Rather, the important result specific to our 
theory is that the divergence of the fluctuation is caused 
by the Goldstone mode for the instanton. We expect that 
this mechanism should be taken into account by combin- 
ing a singular perturbation method with a path integral 
approach for the analysis of microscopic models. 
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